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Abstract 

An Abelian sandpile model is considered on the Husimi lattice of triangles with 
an arbitrary coordination number q. Exact expressions for the distribution of height 
probabilities in the Self-Organized Critical state are derived. 
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1 Introduction. 



To explain the temporal and spatial scaling in dynamical dissipative systems being in the 
Self-Organized Critical (SOC) state, Bak, Tang, and Weisenfeld have introduced cellular 
automaton models known as sandpiles [|TJ. 

Latter on, Dhar has studied sandpile models which can be described by the Abelian 
group 0. For these Abelian sandpile models (ASM's) some ensemble-average quantities 
such as the total number of configurations in the SOC state, distribution of height prob- 
abilities, two-point correlation functions, etc., have been calculated exactly on the square 
lattice [3-6]. 

Unfortunately, for the sandpile models formulated on the standard plane lattices it 
seems to be rather difficult to formulate a direct theoretical approach leading to the 
solution in closed form. 

An appreciable insight can be obtained, however, by examining pseudolattices such as 
the Bethe lattice, Husimi lattice, and their generalizations. 

Using this approach, Dhar and Majumdar J7| have obtained exactly a more complete 
set of sandpile characteristics for the Bethe lattice. 

The aim of this article is to analyze the distribution of height probabilities of the ASM 
on the Husimi lattice with an arbitrary coordination number q. We show that this model 
can be solved exactly by using the fractal structure of the lattice. 

The outline of the article is as follows: in the next section, we define the lattice and 
the ASM on it. In section 3, we consider the Husimi lattice with q = 4 and find recursion 
relations for numbers of allowed configurations in the SOC state. In section 4, we compute 
exactly the distribution of height probabilities in this particular case. In section 5, we 
generalize our results to an arbitrary coordination number of the lattice. 

2 Lattice and Model. 

As a basic building block for constructing a cactus or a pure Husimi tree ||, we will take 
7 triangle plaquettes joined at the base site (root), see Fig. la. This basic building block 
will be called the first- generation branch. To the second- generation branch, we attach 27 
basic building blocks at each free site of the first-generation branch (Fig. lb). Continuing 
this process we develop higher-generation branches. Then, at the final step, we take three 
nth-generation branches and connect their base sites by the triangle plaquette. As a 
result, we will get the cactus with the coordination number q = 2(7 + 1). 

Let as define ASM on this connected graph of iV sites as follows: to each site i 
(1 < i < N) we associate an integer Zj (1 < z% < q) which is the height of a column of 
sand grains. The evolution of the system is specified by two rules: 

i. Addition of a sand grain at a randomly chosen site i increases Zj by 1. 

ii. The site i topples if the height z, exceeds a critical value z c = q and sand grains 
drop on the nearest neighbors. 
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The number of surface sites of the Husimi tree is comparable with the interior ones. 
Hence, the calculation of the thermodynamic limit of the bulk properties requires special 
care. In our work we define the height distribution of sand grains for the sites deep inside 
the tree. Using these interior sites one can construct an infinite lattice, for they have the 
same features. Therefore, we will consider the problem on the Husimi lattice rather than 
on the Husimi tree. 

Any configuration {z{\ on the Husimi tree in which 1 < Zi < q is a stable configuration 
under the toppling rule. These configurations can be divided into two classes: allowed 
and forbidden configurations @. 

In the SOC state, only allowed configurations have a nonzero probability. Any sub- 
configuration of heights F on a finite connected set of sites is forbidden if 

Zi < qi, 

where qi is a coordination number of a site % in the given sub configuration F [0]. 

In turn, we can divide the allowed configurations on an nth-generation branch of the 
Husimi tree into three nonoverlapping classes: weakly allowed of the type 1 (W\), weakly 
allowed of the type 2 (W 2 ) and strongly allowed (S) configurations. 

Consider an allowed sub configuration C on the nth-generation branch G n with a root 
a (Fig. 2a). The coordination number of the root a is q — 2. Adding a vertex b to the G n , 
one defines a subgraph G' = G n U b. If the subconfiguration C = C U b with zt, — 1 on 
the G' is forbidden, then C is called the weakly allowed subconfiguration of type 1 (W\). 
Thus, W\ can be locked by one bond. 

Now add two vertices b and d to the G n and consider a subconfiguration C" = CUbUd 
on G" = G n U b U d (Fig. 2b) . If the subconfiguration C" = CUbUd with z\, = 1, zj = 1 on 
the G" is forbidden, then C is called the weakly allowed subconfiguration of type 2 (W 2 ). 

Any allowed sub configurations defined on the nth-generation branches that cannot be 
locked by one bond or by two bonds form a strongly allowed class (S). 

It is important to note that any subconfiguration of the W\ type is also of the W 2 
type. To obtain the nonoverlapping classes, we always check the subconfiguration first to 
belong to the W\ type and only after that to the W 2 type. 

We will start with mentioning some forbidden subconfigurations that can occur on the 
Husimi tree. 

Let us first consider the subconfiguration shown in Fig. 3a. It is easy to check that 
C = C\ U C 2 will be forbidden if both C\ and C 2 are of the W\ type. The next example 
of the forbidden subconfiguration is shown in Fig. 3b where configurations Ci, C 2 and C3 
should be of the W\ or W 2 type. 
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3 Recursion relations. 

In this section, for simplicity, we will consider a Husimi tree with 7 = 1 (q = 4). The nth- 
generation branch G n with a root a (Fig.4) consists of two (n — l)th-generation branches 
Gr^li and G^-i with roots a\ and respectively. Vertices a\ and a 2 are nearest neighbors 
of the root a. Let Nw x {G n ,i), Nw 2 (G n ,i) and Ns(G n ,i) be the numbers of distinct W\, 
W 2 and S type subconfigurations on the G n with a given height z a = i. Let us also 
introduce 

N Wl (G n )=f2 N w 1 (Gn,i), (3.1) 
i=i 

^(GO =E JV Wi( G! n,i)> (3-2) 



i=i 
4 



iV s (G n ) = £> s (G n ,i). (3.3) 
i=i 

The number of allowed subconfigurations on the G n with a given root is 

N(G n ,i) = N Wl (G n ,i) + N W2 (G n ,i) +N s (G n} i). (3.4) 

At the same time, N Wl (G n ), N W2 (G n ) and N s (G n ) can be expressed in terms of the full 
numbers of allowed subconfigurations on the two (n — l)th-generation branches G^-\ and 

r m . 

^n-l- 

N Wl (G n ) = NsiG^NsiG^) + Ns^G^Nw, (G^) + N Wi (G { 2i)N s (g£1 1 ) + 
NsiG^Nw^l,) + Nw^G^NsiG^) + N Wl (G^_ 1 )N W2 (G^_ 1 )+ (3.5) 

N W2 (G n ) = Ns(G^l 1 )Ns(G^l 1 ) + ^G^AW^-i) + N Wi (g£1 1 )N s (G { *1 1 )+ 

Ns(Gn-i)N w , 2 (Gnl 1 ) + N w . 2 (G ( r ll 1 )N s (Gnl 1 ) + N Wl (Gill 1 )N W2 (G i *l 1 )+ (3.6) 
^w 2 {Gn-i)Nw 1 {G^_i) + N W2 (Gn-i)N W2 (Gn-i), 

N s (G n ) = 2N s (G { n 1) _ 1 )N s (G^ 1 ) + N s (G^l 1 )N Wl (G^l x ) + N Wl (G^ 1 )N s (G ( ^l 1 )+ 

2N s (g£1 1 )N W2 (G { *1 1 ) + 2N W2 (G^1 1 )N S (G { *1 1 ). (3.7) 

As we can see, expressions ( |3.5| ) and ( |3.6| ) are equal to each other. We will prove in 
section 5 that this holds for all q 
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Let us define 

where N w = N Wl = N W2 . 

From (|3.5|) - (|3.7|) one may find the following recursion relation: 

„ lr , l± mG'iU) + 2X(ggQ + 3X(G l ^L 1 )X(G^'_ l ) 
2 + 3X(G< 1 I 1 )+3X(G» 1 ) 



If we consider graphs G^-i and <j^2i to be isomorphic, then N(G^}_ 1 ) = N{G^-i) and 
equation ( |3.9| ) simplifies to 

X(G n ) = ~(l + X(G n _ 1 )) . (3.10) 

With the initial condition X(G ) = 1/2, equation ( |3.10| ) has a simple solution 

X(G n ) = 1 - 2~ (n+1) . (3.11) 

On the infinite generation branches (ri — > oo) the ratio of the number of the Wi type 
configurations or of the type to the strongly allowed ones tends to 1. 



4 Distribution of height probabilities. 

Consider now a randomly chosen site O deep inside the Husimi tree (Fig. 5). The nearest 
neighbors of the site O are roots of the four nth-generation branches. For a given value 
zq = i (1 < i < 4) at the site O, the number of allowed configurations N(i) is expressed 
via iV(G«), N(G^), N(G®), N(G$), where again N(G&) = Et =1 N(G£\i), a = 
1,2,3,4. 

If % — 1, then each allowed configuration on the branches G^ 1 -*, G^ 2 -*, G^\ G^ cannot 
be of the W\ type. It is also evident that two W 2 type configurations cannot occur on the 
neighboring branches Gffl and G^ or G^ and G^\ Hence, the number N(l) of allowed 
configurations with isomorphic branches is 

4 

N(l) = [1 + 4X + 4X 2 ] J] N s (G {a) ). (4.1) 

a=l 

Arguing similarly, one can get the following expressions for the numbers of allowed 
configurations for a given value at the site i 

4 

N(2) = [1 + 8X + 12X 2 ] J] N s (G (a) ) , (4.2) 

a=l 



4 



jV(3) = [1 + 8X + 22X 2 + 12X 3 ] J[ N s (G {a) ) , (4.3) 



Q = l 

4 



JV(4) = [1 + 8X + 22X 2 + 24X 3 ] J[ N s {G {a) ) . (4.4) 



a=l 



The probability P(z) of having the height % at site O is 

-W total 

where X to t a i = Ht=i is the total number of allowed configurations on the Husimi 
tree. 

The system can reach the SOC state only in the thermodynamic limit. For the sites 
far from the surfaces in this limit (n — > oo) we have X = 1. Thus, from (|4.1|)-(|4"75|) we get 



q 21 43 55 

P ( 1)= 128' P < 2) = T28- P(3) =128- P(4) =128' ^ 

The obtained values ( f4.6|) for the height probabilities characterize the SOC state of 
the ASM on the Husimi lattice with q = 4. 



5 Generalization to the Husimi lattice with an arbi- 
trary coordination number q. 

We now intend to generalize results obtained in two previous sections. Before the recursion 
relations and expressions for the height probabilities will be written out, we want to prove 
that the number of the W\ type configurations and the number of the W2 type ones is 
equal on the nth-generation branch of the Husimi tree. 

The numbers of weakly allowed configurations of both the types can be written as 

N Wl (G n )= N Wl (G n ,l)+ N Wl (G n ,2) + .-- + N Wl (G n ,q), (5.1) 
N W2 (G n )= N W2 (G n ,2) + --- + N W2 (G n ,q-l) + N W2 (G n ,q). (5.2) 

It is easy to see that Nw 1 (G n ,l) = Nw 2 (G n ,q) = 0. By definition, the W\ type 
configurations can be locked by one bond. When the height of the root is increased 
by 1, it becomes of W2 type. Similarly, decreasing the height of the root of W2 type 
configuration by 1, we get the W\ type. Thus for each configuration of the W\ type there 
is a unique configuration of the W2 type and vice versa. Therefore, we can conclude 
that (|5.1|) and (|5.2|) are equal. 

For an arbitrary coordination number of the Husimi lattice, equation ([3.10 ) generalizes 

to 

X(G n ) = ^-(l + X(G n - 1 )). (5.3) 
27 



5 



and with the initial condition X(Gq) = 1/(27), it has the solution 



For large enough n we have X(G n ) n =°° 1/(27 ~~ !)■ 

In the same way as for the q = 4 case, we can calculate the height probabilities of the 
sites deep inside the lattice 



P(i) 



(2 7 + 1) 7+1 



27+2( 7 + l)2( 7 + 3) 7 



(i-l)/2 i-odd 
i-l (i-2)/2 i-even / , i\| 

y y t u ; :, — = — -2 ni 3 n2 (2 7 + i)- (rii+ri2) , (5.5) 

where ni + n 2 < 7 + 1 and 1 < i < q. 

The probability of the height to be equal 1 has a simple form 

(27 + 1) 7+1 

P(1) = 27+2( 7+ l)2( 7 + 3 )7 ' (5 ' 6) 

and tends to when 7 — > 00. 

In this article, we have exactly calculated the distribution of height probabilities in the 
Self- Organized Critical state of the ASM on the Husimi lattice of triangles. The next step 
of our investigations will be the Husimi lattice with square plaquettes where we intend 
to derive exact expressions for two-point correlation functions and critical exponents of 
avalanches. The choice of the Husimi lattice gives us hope to find the relationship between 
the chaos and SOC state, as some spin models formulated on this lattice show the chaotic 
behavior H, [TD . 
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(a) 



(b) 




base site base site 



Figure 1: (a) A first-generation branch which consists of 7 triangle plaquettes joined at 
the base site, (b) A second-generation branch. 



(a) (b) 




Figure 2: (a) A nth-generation branch G n and vertex b form a subgraph G'. (b) Now two 
vertexes b and d and the G n form a subgraph G" . 
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(a) (b) 




Figure 3: Two examples of forbidden sub configurations that can occur on the Husimi 
tree. 




a 



Figure 4: The nth-generation branch G n with two nearest (n — l)th-generation branches 
G^U and G^. 
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Figure 5: A site O with height % is located deep inside the lattice and surrounded by the 
four nth-generation branches G^\G^\ G$. 
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